Synthesis and characterization of entangled mesoscopic superpositions 

for a trapped electron 
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We propose a scheme for the generation and reconstruction of entangled states between the 
internal and external (motional) degrees of freedom of a trapped electron. Such states also exhibit 
quantum coherence at a mesoscopic level. 
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A single electron trapped in a Penning trap is one of 
the most fundamental quantum systems. Among its pe- 
culiar features, it allows the measurement of fundamental 
physical constants with striking accuracy. Recently, for 
instance, the electron cyclotron degree of freedom has 
been cooled to its ground state, where the electron may 
stay for hours, and quantum jumps between adjacent 
Fock states have been observed M. It is therefore evi- 
dent that the manipulation and the characterization of 
the state of a trapped electron is an important issue, with 
implications in the very foundations of quantum mechan- 
ics. Earlier works ||pj have dealt with this problem for 
one (motional) degree of freedom. On the other hand, 
entanglement [|| has been recognized as one of the most 
puzzling features of quantum mechanics, being also the 
basis of quantum information processing A striking 
achievement in this rapidly expanding field has been the 
recent entanglement of four trapped ions 0. However, it 
is also possible (and conceptually equivalent) to entangle 
different degrees of freedom of the same particle || . 

In the present work we propose to generate entangled 
states (combined cyclotron and spin states) of an elec- 
tron in a Penning trap by using suitable applied fields. 
The complete structure of the cyclotron-spin quantum 
state is then obtained with the help of a tomographic 
reconstruction from the measured data. 

In a Penning trap an electron is confined by the combi- 
nation of a homogeneous magnetic field along the positive 
z axis and an electrostatic quadrupolc potential in the xy 
plane [[l] . The spatial part of the electronic wave function 
consists of three degrees of freedom, but neglecting the 
slow magnetron motion (whose characteristic frequency 
lies in the kHz region), here we only consider the axial 
and cyclotron motions, which are two harmonic oscilla- 
tors radiating in the MHz and GHz regions, respectively. 
On the other hand, the spin dynamics results from the 
interaction between the magnetic moment of the electron 
and the static magnetic field, so that the free Hamilto- 
nian reads as Q| 



|C0 = Tiuj z a\a z + %uj c a\a c + huj s a z /2 , 



(1) 



Here, in addition to the usual trapping fields, we con- 
sider an external radiation field as a standing wave along 
the z direction and rotating, i.e. circularly polarized, in 
the xy plane with frequency Q ||. To be more specific, 
we consider a standing wave within the cylindrical cavity 
configuration fill with the (dimensionless) wave vector 
k. Then, the interaction Hamiltonian reads j^] 

H iut = he [a c e int + oje - *"*] cos(nz + <f>) 
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where a± = (a x ± ia y )/2, and 



), (2) 
The phase 

ah defines the position of the center of the axial motion 
with respect to the wave. Depending on its value the 
electron can be positioned in any place between a node 
(<f> = 0) and an antinode (ah = ±7r/2). The two cou- 
pling constants e and £ are proportional to the ampli- 
tude of the applied radiation field. For our purposes we 
also consider the possibility to introduce pulsed standing 
waves through the microwave inlet ||] so that e and 
C become time dependent. The duration of the pulse is 
assumed to be much shorter than the characteristic axial 
period, which is of the order of microseconds. Depend- 
ing on f2 and ah, the interaction Hamiltonian (^) gives rise 
to different contributions at leading order in the Taylor 
expansion of sui(kz + ah) and cos(kz + 4>). 

Now, nonclassical cyclotron states can be entangled 
with the spin states through the following steps. First, 
we consider = 0, Q = uj s , and a pulsed standing wave 
of duration Aii = t\ — to = t\. In the following we shall 
work in a frame rotating at the frequency lo s . Then, 
the total Hamiltonian (H[ ICC + H lnt ) can be written as 
H = h,uj z a\a z + h(u) c — u s )a],a c + h((t)a x K(a z + a\). The 
first two terms can be neglected during the pulse dura- 
tion. In fact, the latter is assumed to be much smaller 
than lu^ 1 and (lu c 



is)- 1 » lj- 1 0. Its effect on the 
axial degree of freedom can be described by means of the 
relation p z (t\) 
and p z = —i(a z 



p z (0) - (a x , where C = « dtC(t), 



where the indices z, c, and s refer to the axial, cyclotron 
and spin motions, respectively. 



Subsequently, we allow a free evo- 
lution for a time At 2 = t 2 —t 1 =%/ (2u z ). That amounts 
to having z(t 2 ) = p z (h) = p z (0) - C<r x . 

Finally, we consider the action of another pulsed stand- 
ing wave with ah = — ir/2, Q = uj c , for a time At^ = 



1 



tj, — In such a case the effective Hamiltonian (in 
the frame rotating at the frequency u> s ) becomes H = 

kz. Again, 



Tiuj z a\a z + h(uj c — Lu s )a\a c + he(t)(a c + aft 



since w c — w s is of the order of lj z Q , we can neglect the 
first two terms in the previous Hamiltonian. Then, the 
time evolution operator is equivalent to a displacement 
operator Q U(At 3 ) = D (-iez), with e = k f* 3 dte(t). 
Since z is not affected by the time evolution under the 
previous Hamiltonian, it remains unaltered during Ai 3 , 
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that is U (£3 — to) = cxp j — ie(a, 

means that the state evolution for the spin-cyclotron sys 
tern is given by 



f>(h) = Tr x {£> {-iep z ) D (ieC,o x 

x R(0)ft £>H-i®,)} , (3) 

where R(0) is the initial density operator for the whole 
system. If we consider for instance the initial axial state 
as a Gaussian state with momentum width d, the above 
equation can be rewritten as 



l 



dp z D {-Hp. 



) D (ie(a x ^ 



c/3(0)£t (iKa x ) & {-ilp z ) eM-P 2 Jd 2 ) ■ (4) 



Assuming ( > (!, which is easily obtained in the case 
of the ground state of the axial oscillator, we can ap- 
proximate the evolution of an initial pure state p(0) = 
|$o)($o| hito |$) = D(aa x )\$o), where a = ie(. 

It is now immediate to see that an initial state |$o) = 
1 0)| t) evolves into 



l*> = (lt)|a H 



-»/V2, 



(5) 



where | |) and | J.) denote spin eigenstates, while \a±) = 
M(\a) ± | — a)) are the even-odd coherent states |L3| of 
the cyclotron mode with J\f a normalization factor. Since 
the latter are orthogonal states, Eq. (^|) represents a max- 
imally entangled state H] . On the other hand, even and 
odd coherent states may represent the basis of uncon- 
ventional quantum bits |l5j , hence states of the form of 
Eq. (|^) could be of importance to encoding and manip- 
ulating quantum information Furthermore, a simple 
spin rotation (as it can be seen below) is sufficient to 
realize the transformation 



1$) 



(HIT> + |-a)U»/V2, 



(6) 



a state already discussed in Refs. Jsjjlq]. The electronic 
state (|(]) possesses two very interesting features: first, if 
\a\, i.e. the product e£ is much larger than 1, \ty) is a 
typical example of Schro dinger- cat state Second, 
the full state of the trapped electron is an entangled state 
between the spin and cyclotron degrees of freedom. It is 
worth noting that states like (||) or (||) persist for many 



cycles since the decoherence (i.e. the rapid destruction 
of superposition states due to the entanglement with the 
environment JT^]) in such a system is quite small [Q (in 
contrast, in Ref. fa] the mesoscopic superposition appears 
only cyclically). 

The most general pure state of the trapped electron 
can be cast in the form |^) = ci|^i)| f) + Pa | V2) | I) > 
\ipi) and ^2) being two unknown cyclotron states, and 
the complex coefficients c\ and ci satisfying the normal- 
ization condition | cx | 2 + | C2 1 2 = 1. The density operator 
p = |*)( V I / | associated to the pure state can be ex- 
pressed in the form 



P = 



| Cl | 2 |Vl)«>l| C lC ^l)^2 
C2Ct\^ 2 )(H |C 2 | 2 |V>2><V>2 



(7) 



whose elements p^* — CiC*\ipi)(ipj\ are operators in the 
cyclotron Hilbert space. The phase-space description cor- 
responding to (Q) is given by the Wigner-function ma- 
trix |llS,19i whose elements are 



Wij(a) = (Sij(a - a)) = Tr[/3%(a - a)] , 



(8) 



where 5ij (a — a) is the Fourier transform of the displace- 
ment operator with i,j — 1,2 denoting the spin compo- 
nents. 

In order to characterize the generic state (||) we use 
a simple reconstruction procedure: Adding a particular 
inhomogeneous magnetic field — known as the "magnetic 
bottle" field [Q — to that already present in the trap, it is 
possible to perform a simultaneous measurement of both 
the spin and the cyclotronic excitation numbers. The 
useful interaction Hamiltonian for the measurement pro- 
cess is then (l) 



bottle 



TlUJh 



it; 



9 o z 



a'a r H 

c 2 2 



(9) 



where the angular frequency ujb is directly related to the 
strength of the magnetic bottle field. 

Eq. (^J) describes the fact that the axial angular fre- 
quency is affected both by the number of cyclotron ex- 
citations n c = a\d c and by the eigenvalue of a z . The 
modified (shifted) axial frequency can be experimentally 
measured |l| after the application of the inhomogeneous 
magnetic bottle field. One immediately sees that it as- 
sumes different values for every pair of eigenvalues of 
h c and a z , due to the fact that the electron g factor 
is slightly (but measurably [[I)) different from 2. Then, 
repeated measurements of this type allow us to recover 
the probability amplitudes associated to the two possible 
spin states and the cyclotron probability distribution in 
the Fock basis. The reconstruction of the density matri- 
ces \il>i)(i(>i\ (i = 1, 2) in the Fock basis is then possible by 
employing a technique similar to the Photon Number To- 
mography (PNT) which exploits a phase-sensitive 
reference field that displaces in the phase space the par- 
ticular state one wants to reconstruct 



2 



Following Ref. [0, immediately before the measure- 
ment, we apply a pulsed standing wave (^) tuned to 
Q = w c with <p = Q in order to get a displacement 
7 = —ie/n on the cyclotron. Thus we can interpret the 
quantity 



pW(n,7) = (n|L» t ( 7 )p^£i(7)|n) 



(n, 7 |p("V,7) 



(10) 



as the probability of finding the cyclotron state \ipi) in 
a displaced number state |n, 7) [p2| . It should be noted 
at this stage that the probability distribution (|l0|) is not 
normalized to unity. Instead, one has 



00 

E 

n=0 



P W (n,7) = N 2 



(11) 



Eq. ([ll]) allows us to retrieve in a simple way the moduli 
of the coefficients c\ and C2 in Eq. (0) from the measured 
data. Fixing a particular value of 7, it is then possible 
to recover the probability distribution ( |l0|) performing 
many identical experiments. 

Expanding the density operator p( u > in the Fock basis, 
and defining 7V C as an appropriate estimate of the maxi- 
mum number of cyclotronic excitations (cut-off), we have 



N, 

E 

fc,m— 



P«(n, 7 )= E (n,-y\k)(k\fiM\m)(m\n,>y) . (12) 



The projection of the displaced number state |n, 7) onto 
the Fock state \m) can be obtained generalizing the result 
derived in Ref. |23[| . 

Let us now consider, for a given value of |7|, pW(n, 7) 
as a function of ip = arg[ 7 ] |2j| and calculate the coeffi- 
cients of the Fourier expansion 



pW(n,s) 
for s = 0,1,2,.. 

>(0 



27T ,/ 



Combining Eqs. (ph and (|R 



7V c -s 



pW(n, S )= J] Gi s ) m (|7|)(m + ,s|^)|m) 



m=0 



(13) 
we get 

(14) 



where the explicit expression of the matrices G is given 
in Ref. g|. 

We may now note that if the distribution (n, 7) 
is measured for n S [0, AT] with AT > iV c , then Eq. Jl4| ) 
represents for each value of s a system of Af + 1 linear 
equations between the A^+l measured quantities and the 
N c + 1 — s unknown density matrix elements. Therefore, 
in order to obtain the latter, we only need to invert the 
system 



s\p {ii) \m) 



N 

E 

n=0 



M« (| 7 |)P«(tm), (15) 



where the matrices M are given by M = (G 1 G) 1 G 



Since the overdetermined system ( |l4|) is inverted using 
the method of least squares, we are sure that when the 
measured probabilities are slightly inaccurate, the quan- 
tities calculated from the reconstructed density matrix 
best fit the measured ones 

However, this kind of measurement does not allow to 
retrieve the relative phase between the complex coef- 
ficients ci and C2- We can then again use the Hamilto- 
nian (g) tuned to have a spin rotation, i.e. f2 = ui s and 
4> = tt/2. After a n/2 spin rotation we have 



|*)-^|*) = ^[( Cl |^>-iP2|^))|T) 
+Hd|Vl>+C2|V2»| i>] 



(16) 



We can now repeat the spin measurements just as we 
have described above in the case of the unknown initial 
state I "J). Repeating this procedure over and over again 
(with the same unknown initial state) for a large number 
of times and tracing out the cyclotron degree of freedom 
(no drive in this case is required), it is possible to recover 
the probabilities P^ associated to the two spin eigenval- 
ues for the state \^). Without loss of generality, we can 
assume ci S R, C2 = |c2|e l9 , and (V'ilV'2) = re 4 ' 3 , which 
yield 



2r| Cl ||c 2 |sm(0 + /3)] 



P(D = J[1 



P {2) = i[l-2r|ci||c 2 |sin(0 + /3)] . 



(17a) 
(17b) 



It is important to note that the probabilities P^ can 
be experimentally sampled and that the modulus r and 
the phase of the scalar product (tpilfa) can be both 
derived from the reconstruction of the cyclotron density 
matrices o^ 11 ) and p( 22 ) It 



Thus we are able to find 



the relative phase 9 by simply inverting one of the two 
Eqs. @, e.g. 



2P( J 



1 



2r ci c 2 



■0 



(18) 



where the ambiguity in the arcsin function can be re- 
moved by repeating the procedure above using a second 
spin rotation ( |l6[ ) with an angle different from tt/2. As in 
any tomographic scheme, our reconstruction procedure 
yields the state to be measured up to an uninteresting 
phase factor which, however, does not affect the results. 

As an example of the states that can be generated with 
the method outlined above, and of the application of the 
proposed reconstruction procedure, we show in Fig. [l]the 
results of numerical Monte-Carlo simulations of the re- 
construction of the state @ . In this simulation we have 
used the value \a\ = 1.5 which is experimentally accessi- 
ble [jlj and gives mesoscopic entangled superpositions of 
cyclotron coherent states with opposite phase. In order 
to account for actual experimental conditions, we have 



3 



also considered the effects of a non-unit quantum effi- 
ciency rj in the counting of cyclotron excitations. When 
r\ < 1, the actually measured distribution is related to 
the ideal distribution by a binomial convolution p(| . As 
it can be seen from Fig. [I], the reconstructed distribu- 
tions turn out to be quite faithful. We also would like 
to emphasize that the particular shape of W12 is due to 
the quantum interference given by the entanglement be- 
tween the two degrees of freedom: in fact, in absence of 
entanglement p^ 12 ^ would just be a replica of the diagonal 
parts p( n ) and p( 22 K 



True Reconstructed 




FIG. 1. Simulated tomographic reconstruction of the 
Wigner matrix Wij = Wij/c t c* for the state of Eq. (^) with 
a — 3i/2. The quantum efficiency is r\ — 0.9 and 10 6 data per 
phase have been simulated. In this simulation an amplitude 
I7I = 1.2 of the applied reference field has been used. The 
theoretical distributions are plotted on the left for compari- 
son. 

We have performed a large number of simulations with 
different states and several values of the parameters, 
which confirm that the present method is quite stable 
and accurate. In addition, and for all the cases consid- 
ered, the values of the parameters cy, C2, and 9 are very 
well recovered, with a relative error of the order of 10~ 5 . 

To conclude, we have proposed here a method which 



is able to synthesize and characterize highly nonclassical, 
maximally entangled states for a single trapped electron. 
The method is based on simple operations (switching on 
and off standing waves) which are currently realized in 
the present trap technology An experimental imple- 
mentation of the proposed method might yield new in- 
sight in the foundations of quantum mechanics and allow 
further progress in the field of quantum information. 
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